3D Relativistic Hartree-Bogoliubov model with a separable 

pairing interaction 

T. Niksic and D. Vretenar 
Physics Department, Faculty of Science, 
University of Zagreb, 10000 Zagreb, Croatia 

Yuan Tian and Zhong-yu Ma 
China Institute of Atomic Energy, Beijing 102413, People's Republic of China 

P. Ring 

Physik- Department der Technischen Universitdt Miinchen, D-85748 Garching, Germany 

(Dated: December 7, 2009) 

Abstract 

A recently introduced separable pairing force for relativistic Hartree-Bogoliubov (RHB) calcu- 
lations, adjusted in nuclear matter to the pairing gap of the Gogny force, is employed in the 3D 
RHB model for triaxial shapes. The pairing force is separable in momentum space but, when 
transformed to coordinate space in calculations of finite nuclei, it is no longer separable because of 
translational invariance. The corresponding pairing matrix elements are represented as a sum of 
a finite number of separable terms in the basis of a 3D harmonic oscillator. The 3D RHB model 
is applied to the calculation of binding energy surfaces and pairing energy maps for a sequence of 
even-A Sm isotopes. 

PACS numbers: 21.60.Jz, 21.60.Ev, 21.10.Rc, 21.10.Ky 
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I. INTRODUCTION 



The structure of heavy complex nuclei with a large number of active valence nucleons is, 
at present, best described by the framework of nuclear energy density functionals (NEDF). 
A variety of structure phenomena, not only in stable nuclei, but also in regions of nuclei 
far from the valley of /^-stability and close to the nucleon drip-lines, have been described 
with self-consistent mean-field models based on the Gogny interaction, the Skyrme energy 
functional, and the relativistic meson-exchange effective Lagrangians Q| • 

Self-consistent relativistic mean-field (RMF) models have been employed in analyses of 
properties of ground and excited states in spherical and deformed nuclei. For a quantitative 
analysis of open-shell nuclei it is necessary to consider also pairing correlations. Pairing 
has often been taken into account in a very phenomenological way in the BCS model with 
the monopole pairing force, adjusted to the experimental odd-even mass differences. In 
many cases, however, this approach presents only a poor approximation. The physics of 
weakly-bound nuclei, in particular, necessitates a unified and self-consistent treatment of 
mean-field and pairing correlations. This has led to the formulation and development of 
the relativistic Hartree-Bogoliubov (RHB) model, which represents a relativistic extension 
of the conventional Hartree-Fock-Bogoliubov framework. In most applications of the RHB 
model the pairing part of the Gogny force has be employed in the particle-particle 
(pp) channel: 



j=l,2 



with the set DIS 



^ for the parameters Wi, Bi, Hi, and Mj (i = 1, 2). A basic advantage 
of the Gogny force is the finite range, which automatically guarantees a proper cut-off 
in momentum space. However, the resulting pairing field is non-local and the solution 
of the corresponding Dirac-Hartree-Bogoliubov integro-differential equations can be time- 
consuming, especially in the case of deformed nuclei. An alternative is the use of a zero- 

n 

range, possibly density-dependent, (5-force in the pp-channel of the RHB model [5|, but 
this approach introduces an additional cut-off parameter in energy. The effective range 
of the pairing interaction is determined by the energy cut-off, and the strength parameter 



must be chosen accordingly in order to reproduce empirical pairing gaps. In Ref. 



6j we have 



implemented a renormalization scheme for the relativistic Hartree-Bogoliubov equations with 
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a zero-range pairing interaction. The procedure is equivalent to a simple energy cut-off with 
a position dependent coupling constant, and the resulting average pairing gaps and pairing 
energies do not depend on the cut-off energy. A density-dependent strength parameter of 
the zero-range pairing can be adjusted in such a way that the renormalization procedure 
reproduces in symmetric nuclear matter the pairing gap of the Gogny force Eq. ([T]). 

In a series of recent articles we have introduced a separable form of the pairing 

force for RHB calculations in spherical and axially deformed nuclei. The force is separable 
in momentum space, and is completely determined by two parameters that are adjusted 
to reproduce in symmetric nuclear matter the bell-shape curve of the pairing gap of the 
Gogny force. Technically this approach is similar to the derivation of non-empirical pairing 



functionals from low-momentum interactions 
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In applications to finite nuclei, when 



transformed from momentum to coordinate space, the pairing force is no longer separable 
because of translational invariance. It has been shown, however, that a method developed 
by Talmi and Moshinsky can be used to represent the corresponding pp matrix elements as 
a sum of a finite number of separable terms. When the nucleon wave functions are expanded 
in a harmonic oscillator basis, spherical or axially deformed, the sum converges relatively 
quickly, i.e. a reasonably small number of separable terms reproduces with high accuracy 
the results of calculations performed in a complete basis. 

The simple separable force considered in Refs. 0, S [s] reproduces pairing properties of 
spherical and axially deformed nuclei calculated with the original Gogny force, but with the 
important advantage that the computational cost is greatly reduced. In the present work we 
extend this approach to triaxial nuclei, and introduce a 3D RHB model with the separable 
pairing interaction of Ref. j3] in the pp channel. This model will enable systematic calcula- 
tions of RHB binding energy surfaces in the /3 — 7 plane, based on relativistic nuclear energy 
density functionals, that can be used as input for the generator coordinate method con- 
figuration mixing of angular-momentum projected triaxial wave functions, or to determine 
the parameters of a five-dimensional collective Hamiltonian for quadrupole vibrational and 
rotational degrees of freedom 12|. In Sec. [Ill we introduce the 3D RHB model and derive 
the pp matrix element of the pairing force as a sum of a finite number of separable terms in 
the basis of a 3D harmonic oscillator. Illustrative calculations for even-A Sm isotopes are 
presented in In Sec. IIIII Section IIVI summarizes the results and ends with an outlook for 
future applications. Details on the expansion of single-nucleon spinors in the 3D harmonic 
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oscillator basis, and the transformation of the product of harmonic oscillator wave functions 
to the center-of-mass frame, are included in the two appendixes. 



II. 3D RELATIVISTIC HARTREE-BOGOLIUBOV MODEL WITH A SEPARA- 
BLE PAIRING INTERACTION 



The relativistic Hartree-Bogoliubov framework provides a unified description of 
particle-hole (ph) and particle-particle (pp) correlations on a mean-field level by using two 
average potentials: the self-consistent mean field that encloses all the long range ph cor- 
relations, and a pairing field A which sums up the pp-correlations. The ground state of 
a nucleus is described by a generalized Slater determinant |$) that represents the vacuum 
with respect to independent quasiparticles. The quasiparticle operators are defined by the 
unitary Bogoliubov transformation of the single-nucleon creation and annihilation operators: 

(^t = Y.^ik4 + yikCi, (2) 

where U and V are the Hartree-Bogoliubov wave functions determined by the solution of 
the RHB equation. In coordinate representation: 











)l 








[ Vkir) I 




[ Vkir) j 



(3) 



hi) — m — X A 

-A* -h}j + m + A 

In the relativistic case the self-consistent mean-field corresponds to the single-nucleon Dirac 
Hamiltonian h^. In the usual a, u, and p meson-exchange representation, and for the sta- 
tionary case with time-reversal symmetry, i.e. for the ground-state of an even-even nucleus: 

ho = -iaV + /?(m + g^crir)) + g^uj\r) + gpT3p^{r) + e^^A°(r) . (4) 

The classical meson fields are solutions of the stationary Klein-Gordon equations: 

[-A + ml]a{r) = -g^{p)p,{r), (5) 

[-A + ml]u;\r)=gMpir), (6) 

[-A + ml]p\r)=g,{p)pl{r), (7) 

-AA%r) = p,{r) , (8) 

for the a meson, the time-like components of the uj meson and p meson, and the Poisson 
equation for the vector potential, respectively. In the general case when the meson-nucleon 



couplings Qa, guj, and gp explicitly depend on the nucleon (vector) density/?, there is an 
additional contribution to the nucleon self-energy - the rearrangement term , essential for 
the energy-momentum conservation and the thermodynamical consistency of the model. 

In Eq. m is the nucleon mass, and the chemical potential A is determined by the 
particle number subsidiary condition in order that the expectation value of the particle 
number operator in the ground state equals the number of nucleons. The pairing field A 
reads 

Aab{r, r') = ^Yl ^-bcd{r, r')^i,d{r, r'). (9) 

c,d 

where Vabcdi'^^, f') are the matrix elements of the two-body pairing interaction, and the indices 
a, 6, c and d denote the quantum numbers that specify the Dirac indices of the spinor. 
The column vectors denote the quasiparticle wave functions, and Ek are the quasiparticle 
energies. The dimension of the RHB matrix equation is two times the dimension of the 
corresponding Dirac equation. For each eigenvector (f/^, Vk) with positive quasiparticle 
energy Ek > 0, there exists an eigenvector iV^,Ul) with quasiparticle energy —Ek- Since 
the baryon quasiparticle operators satisfy fermion commutation relations, the levels E^ and 
—Ek cannot be occupied simultaneously. For the solution that corresponds to a ground state 
of a nucleus with even particle number, one usually chooses the eigenvectors with positive 
eigenvalues Ek- 

The single-particle density and the pairing tensor, constructed from the quasi-particle 
wave functions 

pUr,r') = J2ycUr)V,k{r'), (10) 

fe>0 

'^cdir,r') = J2u:kir)V,kir'), (11) 

fe>0 

are calculated in the no-sea approximation (denoted by A; > 0): the summation runs over 
all quasiparticle states k with positive quasiparticle energies E^ > 0, but omits states that 
originate from the Dirac sea. The latter are characterized by a quasiparticle energy larger 
than the Dirac gap (^ 1200 MeV). 

Pairing correlations in nuclei are restricted to an energy window of a few MeV around 
the Fermi level, and their scale is well separated from the scale of binding energies, which 
are in the range of several hundred to thousand MeV. There is no empirical evidence for any 
relativistic effect in the nuclear pairing field A and, therefore, a hybrid RHB model with a 
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non-relativistic pairing interaction can be formulated. For a general two-body interaction, 
the matrix elements of the relativistic pairing field read 

(12) 

where the indices (pi,P2,P3,P4: = f,g) refer to the large and small components of the quasi- 
particle Dirac spinors: 

\Wui^.s,t) J yigv{r,s,t) J 

In practical applications of the RHB model to finite open-shell nuclei, only the large com- 
ponents of the spinors f/fc(r) and Vfc(r) are used to build the non-relativistic pairing tensor 
k in Eq. (ITT]) . The resulting pairing field reads 

(14) 



The other components: A/^, A^/, and A^^ can be safely omitted |13 |. 

The Dirac-Hartree-Bogoliubov equations and the equations for the meson fields are solved 
by expanding the nucleon spinors f/(r, s, t) and V{y, s, t), and the meson fields, in the basis 
of a three-dimensional harmonic oscillator (HO) in Cartesian coordinates. In this way both 
axial and triaxial nuclear shapes can be described. In addition, to reduce the computational 
task, it is assumed that the total densities are symmetric under reflections with respect to all 
three planes xy, xz and yz. When combined with time-reversal invariance, this also implies 
that parity is conserved. The single-nucleon basis is defined in Appendix Rl 

In Ref. [3] a new separable form of the pairing interaction has been introduced, with 
parameters adjusted to reproduce the pairing properties of the Gogny force in nuclear matter. 
The gap equation in the ^Sq channel reads 

and the pairing force is separable in momentum space: 

{k\v'^''\k') = -Gp{k)p{k') . (16) 

By assuming a simple Gaussian ansatz p{k) = e"*^^^^, the two parameters G and a have been 
adjusted to reproduce the density dependence of the gap at the Fermi surface, calculated with 



a Gogny force. For the DIS parameterization of the Gogny force: G = —728 MeVfm^ 
and a = 0.644 fm. When the pairing force Eq. (fT6l) is transformed from momentum to 
coordinate space, it takes the form: 

V{r^, r2, r[, r',) =G5{R- R') P{r)P{r')^ (1 - P^) , (17) 

where = | (ri + and r = ri — r2 denote the center-of-mass and the relative coordi- 
nates, and P{r) is the Fourier transform oi p{k): 

P{r) = ^e-^^/^'^^ . (18) 

The pairing force has finite range and, because of the presence of the factor 5 [R — R'), it 
preserves translational invariance. Even though 6 {R — R') imphes that this force is not 
completely separable in coordinate space, the corresponding pp matrix elements can be 
represented as a sum of a finite number of separable terms in the basis of a 3D harmonic 
oscillator: 

Af« A^^ NO 



Kl,. = T.T.Y. ' (19) 

where N^-, Ny, and A^^ are the quantum numbers of the corresponding ID HOs in the center- 
of-mass frame (cf. Appendix [B]) . This means that the pairing field can also be written as a 
sum of a finite number of separable terms: 

AfO atO 

^"/^ = ^ E E E p-^-.^^yS"'"^ > (20) 

N.j,=QNy=QN^=Q 

with the coefficients 

Pn^NyN^ = V^^s''^^^'' i^^s ■ (21) 

The advantage of using the separable pairing interaction Eq. (IT7|) is that the matrices 
yN^NyN:, calculated only once at the beginning of a self-consistent calculation. The 
coefficients PN^NyN^ are re-calculated at each iteration step, using the corresponding up- 
dated pairing tensor k. 

In the following we calculate the antisymmetric matrix element of the pairing interaction 
Eq. (dZD 

(a7| V \P6)^ = (a7| V \P6) - {a^\ V \6P) , (22) 
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in the basis of a 3D harmonic oscillator, where \a) and |a) denote the positive and negative 
x-simplex operator eigenstates (cf. Appendix 



\a) ^ ^ = +> = 1^") ^ 



(23) 
(24) 



The matrix element can be separated into a product of spin and coordinate space factors 



{a^\V\p5)^ = {a^\W-{l-Pn\P5)^ . 
The operator ^(1 — P'^) projects onto the 5 = spin-singlet product state 



(25) 



1 .„S,„I3 



5=0 



= -/^^"-^ (-l)"^""' 1 + {-l)<^< [|Ti) - lit)] |nV> , (26) 



and the problem is reduced to the calculation of the spatial part of the matrix element 



(a7| V \p5)^ = ±^"?+«^+<+<(_i)"?+"S ^1 + (_i) 



1 + (-1)< 



X 



(27) 
(28) 



For W{ri, ra, r[, r'^) = G6{R- R') P{r)P{r') (cf. Eq. ([HD ), the spatial part of the matrix 
element 



0nJ?'l)0n^(?'2)t^(^^l, ^-2, K> ^2)0n/3(K)0n*(^2)'^^l'^^2C?r-W^5 
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can be decomposed into three Cartesian components 



/ = GIJJ, 



(29) 



(30) 



Here we only derive a detailed expression for the x-component 

h = j (pn^{xi,b^)(pnz{x2,b^)P{x)5{X -X')P{x')(p^p{x[,b^)(pni{x'2,b^^ ■ (31) 

We will use the transformation of a product of ID HO wave functions to the center-of-mass 
and relative coordinates, derived in Eq. ( 1B20I1 



^n^{Xl,b,)(f)^.{x2,b,) = J2 Klni't>NSX,B)(t)nSx,b,^r), 
K^i^'l^^^)^ni{x'2,bcc) = ^ M'^^^;(j)N^{X',B)(j)n'^{x',b^,r) 



(32) 
(33) 



The volume element in Eq. (13T!) is transformed: 



dxidx2dx'idx2 = —dXdX'dxdx'. (34) 
The integral Eq. (I3T1) now reads 

4 = - 5^ M:sJ.^ / Pix)<l>^.ix,K,r)dx I Pix')<l>^,ix',K,r)dx'x 

X jj^i^- X')<PnAX, B)<Pn>{X', B)dXdX'. (35) 
By making use of: 

jj KX- X')<PnAX, B)4>n'SX\ B)dXdX' = j <PnAX, B)4>n'SX, B)dX = 6n^,n^ , (36) 
the integral can be written in the form 

L = -J2 M-5=M<„^^/„, {b^,r)In'^ {b,,r) , (37) 

where denotes 

InAbx,r) = J P{x)(j)nSx,h^^r)dx . (38) 

Note that the conditions 

nl + rf^ = n^ + N, and n| + = <. + A^^ , (39) 

have been used to eliminate the sums over ra^ and n'^. To evaluate the integral In^{bx^r), we 
make use of the generating function for the HO wave functions Eq. (lB2p . and calculate the 
following integral: 

JM = j ji-^P^b)Pix)dx = n~^^^^^^^Y.i-ir^^^^ (40) 

Using the definition of the generating function Eq. (IB2p . this expression can also be written 
as follows 

°° / 92" r°° 

J{p,h) = gp2n^__ j p(3;)0^„(3;,5) . (41) 

The series contains only even powers because P{x) is a symmetric function. By comparing 
Eqs. dlO]) and (jH]), we obtain 

£ P(.),„(.» ^ .-V.yi^(_i)»/. (^)"^„,_ . (42) 
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and finally, inserting the relative oscillator length br = V2b (cf. Appendix [B]) , 

^nM,r) = T^^l^^A-ir' <e.en • (43) 



The number of terms in Eq. (1371) and then, of course, in Eqs. (fT9|) and (120|) is in principle 
limited by the dimension of the oscillator basis. If single-particle oscillator states \nxnynz) 
with Tlx + Uy + Uz < N^'"'^ form the basis, the summation over the quantum number of the 
ID HO in the center-of-mass frame in Eq. ( 1371) runs over A'a; = 0, . . . , 2NJ^°'^. This means 
that the maximal total number of terms in Eqs. ( |T9l) and (l20l) equals N"^^^ = (2NJ^"-^ + 1)^. 
However, results of calculations performed in Refs. suggest that the actual number of 

terms that give significant contributions to the pairing field is much smaller. If a cut-off 
condition is imposed 

N.<K, Ny<Nl, Nz<N:, (44) 
the total number of separable terms becomes: 

A^sep = liK + 2) (iV; + 2) {N: + 2) . (45) 

In the next section we will compare some results of illustrative 3D RHB calculations with 
those obtained assuming axial symmetry 9|]. For a meaningful comparison with results 
calculated using the axial RHBZ code 9|], we will make the following choice: A^^^ + Ny < 
and Nz < N^. In this case the total number of separable terms equals 

III. ILLUSTRATIVE CALCULATIONS 

In this section we present the results of illustrative 3D RHB calculations for a sequence 
of Sm isotopes. The separable pairing force Eq. (fT7|l is used in the pp channel, and the 
mean field is determined by the density- dependent meson-exchange effective interaction 
DD-ME2 1^ in the ph channel. DD-ME2 has been adjusted to empirical properties of 
symmetric and asymmetric nuclear matter, binding energies, charge radii, and neutron radii 
of spherical nuclei. The interaction has been tested in calculations of ground state properties 
of large set of spherical and deformed nuclei. An excellent agreement with data has been 
obtained for binding energies, charge isotope shifts, and quadrupole deformations. When 
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used in the relativistic RPA, DD-ME2 reproduces with high accuracy data on isoscalar and 



isovector collective excitations 



14 



In Figs. [T] and [2] we display the self-consistent RHB triaxial quadrupole binding energy 
maps of the i34-i56gj^ isotopes in the /3 — 7 plane (0° < 7 < 60°). The map of the energy 
surface as a function of the quadrupole deformation is obtained by imposing constraints on 
the axial and triaxial quadrupole moments. The method of quadratic constraint uses an 
unrestricted variation of the function 

$^C2,((g2,)-g2,)', (47) 

/i=0,2 

where {H) is the total energy, and (Q2/i) denotes the expectation value of the mass 
quadrupole operators: 

Q20 = - - ?/2 and Q22 = x^-y\ (48) 
q2^ is the constrained value of the multipole moment, and the corresponding stiffness 



constant 



16|. 



The energy maps shown in Figs. [T] and [2] nicely illustrate the gradual transition from the 
prolate and 7-soft deformed light isotopes ^^^■^'^^Sm, through the spherical = 82 neutron 
closed-shell nucleus ^^^Sm, to the strongly prolate deformed, axial nuclei ^^^'^^^Sm. The 
isotopes below the = 82 closed shell are all 7-soft and, just before the shell closure, 
one finds a slightly oblate minimum in ^^°Sm. Sm nuclei with > 82 quickly develop a 
pronounced prolate deformation, much stiffer with respect to the 7 degree of freedom than 
isotopes below A^ = 82. Heavy Sm isotopes are characterized by axially symmetric shapes 
with pronounced prolate minima at /5 > 0.3. 

The binding energy maps correspond to self-consistent solutions of the RHB equations, 
obtained by expanding the nucleon spinors and the meson fields in the basis of a three- 
dimensional harmonic oscillator (HO) in Cartesian coordinates. In the present calculation 
the basis includes N'p"-^ = 14 major oscillator shells. In Figs. [3] and H] we plot the corre- 
sponding contour maps of the proton and neutron pairing energies in the /3 — 7 plane for the 
three lighter isotopes I34,i36,i38gj^^ ^^^^ j^j, ^.j^g three heavier nuclei I52,i54,i56gj^^ respectively. 
Using the separable pairing force Eq. ( fTTj) . the pairing field Eq. ( |20l) is calculated as a sum 
of a finite number of separable terms. 3D calculations in the pp channel have been checked 
by comparing the results for the ground state properties of ^^^~^^^Sm, with those obtained 
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using the axial RHBZ code with the same separable pairing force jol, and with the pairing 
part of the original Gogny force ^, • In the case when axial symmetry is assumed, the 
expansion for the pairing field runs over the quantum numbers A'^^ and Np of the HO in the 
center-of-mass frame, corresponding to the z and p coordinates of the cylindrical coordinate 
system: 



^12 



Gl^l^W^r^P^P^^. (49) 



The maximal values for the quantum numbers in the expansion of Dirac spinors are ?7,° = 
jymaa; ^0 _ jsfj-"-^ i.e. the maximal values for the coefficients in expansion fl49p are 
ATO = 2nl = 2NJ"''^ and = 2ra° = Np"^ . In Ref. ^\ it has been shown that, for axial 
calculations of prolate deformed nuclei, sufficient accuracy is achieved if the expansion of 
pairing matrix elements is limited to: Np < Np = 5, and Nz < N^ = 14. For this choice 
of the cut-off in the expansion of the pairing matrix elements in the basis of the HO in the 
center-of-mass frame, the resulting pairing energies reproduce to a very good approximation 
results obtained with the calculation in the full basis, and also those obtained with the 
Gogny force DIS in the pairing channel. In the present 3D calculation we have, therefore, 
imposed the following cut-off condition for the expansion in Eqs. (fT9l) . (!20|) and (!2T]) : 

iV, = 0, 2, . . . , iV° = 14 and N^ + Ny = 0,2, . . . , 2N^ = 5. (50) 

In Fig. [5] we display the 3D RHB ground-state binding energies for the Sm isotopes 



(72 < N < 92), in comparison with data from the compilation of Audi and Wapstra 18|. 
Calculations have also been performed with the axial RHBZ code j^l, and the inset we plot 
the relative differences (in percent): (^'R-Hbz _^3Drhb^^^3Drhb^ between the corresponding 
ground-state binding energies. As a further test. Fig. E] compares the 3D RHB and axial 
(RHBZ) results for the self-consistent ground-state quadrupole deformations, neutron and 
proton pairing energies of even-A Sm isotopes. In calculations with axial symmetry (RHBZ) 
both the separable force and the Gogny DIS force 4] are used in the pairing channel. The 
excellent agreement between the three sets of results demonstrate not only the numerical 
accuracy of the new 3D computer code, but also show that using the separable pairing force 
in deformed nuclei, virtually identical pairing energies are calculated as with the original 
Gogny force. 
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IV. SUMMARY AND OUTLOOK 



Realistic self-consistent mean-field calculations based on finite-range interactions, includ- 



ing exchange 



challenge 
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;erms and/or pairing correlations, still present a considerable computational 



20 



21 



221 ]. particularly if one considers complex triaxial shapes or extensions 



beyond the simplest mean-field approximation. A great advantage of mean-field models 
based on Skyrme -like zero-range interactions is that they provide a simple and elegant 



treatment of Fock exchange and pairing terms |23|, |2J]. The disadvantage of such forces, i.e. 
the fact that they are constant in momentum space and can induce scattering of nucleons 
very high up into the continuum, does not appear at the Hartree or Hartree-Fock level, at 
which one considers only momenta up to the Fermi surface. In the pairing channel, however, 
because of the specific form of the BCS or Bogoliubov ansatz that takes into account pairing 
correlations on the mean-field level, ultraviolet divergencies occur for zero-range forces. One 
possible solution present the various cutoff procedures that have been used in the literature 
(cf. Ref. 25|] and references therein). These approximations all include additional non- 
physical cut-off parameters. This does not cause any problem in investigations along the 
valley of beta-stability, where gap parameters can be deduced from experimental masses. 
However, the use of cut-off parameters limits the predictive power of such models in un- 
known regions of the nuclear chart, such as for superheavy elements or very neutron rich 
isotopes. 

A completely different approach to the treatment of pairing correlations is the use of 



separable forces. A separable 
has recently been introduced 



brm of the pairing force for RHB calculations in finite nuclei 







The force is separable in momentum space, and is com- 
pletely determined by two parameters that are adjusted to reproduce in symmetric nuclear 
matter the bell-shape curve of the pairing gap of the Gogny force. Because of transla- 
tional invariance, the pairing force is no longer exactly separable in coordinate space, but 
Talmi-Moshinsky techniques allow a simple transformation into a quickly converging series 
of separable terms in a harmonic oscillator basis. Although different from the Gogny force, 
the corresponding effective pairing interaction has been shown to reproduce with high ac- 
curacy pairing gaps and energies calculated with the original Gogny force, both in spherical 
and axially deformed nuclei. In particular, this approach retains the basic advantage of the 
finite-range Gogny force, i.e. the natural cut-off in momentum space. 
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Applications 
deformed nuclei 



lave so far been restricted to the description of spherical [3] and axially 
9| . In this work we have extended the model to describe triaxially deformed 
nuclei. The numerical accuracy of the new model has been analyzed by comparing results 
with those obtained in axially-symmetric calculations, using both the separable force as well 
as the original Gogny DIS force in the pairing channel. 

To illustrate the applicability of this force in the description of realistic triaxial systems, 
we have explored the chain of Sm isotopes with Z = 62 protons, ranging from ^^^Sm to ^^^Sm. 
For the magic neutron number = 82, i.e. for the ^^^Sm isotope, a stable spherical minimum 
is found in the /3 — 7 plane. The two neighboring nuclei ^^^Sm and ^^^Sm are still spherical, 
but with much softer energy surfaces. In heavier isotopes we find a rather rapid transition to 
prolate shapes with a well pronounced minima, and increasing /5-deformation up to /3 ~ 0.3. 
In these heavier nuclei we also find a soft saddle point on the oblate side, that eventually 
becomes a shallow second minimum in the isotope ^^^Sm. Decreasing the neutron number 
below the closed shell at = 82, a 7-soft valley develops with increasing /3-deformation. 
For ^^°Sm isotope we find a shallow oblate minimum, whereas for the lighter isotopes the 
minima are located on the prolate side, but the calculation predicts large fluctuations in the 
7-direction. 

One can envisage many possible applications of the separable pairing force. The force is 
simple enough to be applied in otherwise time consuming calculations, e.g. description of 
triaxial effects, rotating nuclei, fission process, spherical and deformed QRPA. It can also 
be used in various beyond mean-field extensions, such as restoration of broken symmetries, 
fluctuations of quadrupole moment and particle-vibration coupling. In the current version 
of the model, the pairing force has been adjusted to the pairing gap of the phenomenological 
Gogny DIS force. In the next step one could adjust the effective force to a pairing gap 
in nuclear matter calculated in a fully microscopic approach starting from inter-nucleon 



interactions 



26 
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APPENDIX A: THE SINGLE-NUCLEON BASIS 

The Dirac single-nucleon spinors are expanded in the basis of eigenfunctions of a three- 
dimensional harmonic oscillator (HO) in Cartesian coordinates. In one dimension: 

Ki^,) = b;'^'K,H^,{Qe-^'/' , (/i = X, y, z) (Al) 

= and the oscillator length is defined as 



h, = ^j^. (A2) 
" ma;,, 



The normalization factor reads 

Ar„ = 7r"^/^(2"n!)-'/', (A3) 



and Hn{C) denotes the Hermite polynomials [29j 

/oo 
-oo 



(A4) 



The basis state can be defined as the product of three HO wave functions (one for each 
dimension) and the spin factor: 

0^(r;m,) = 0n,(^x)0n,(^y)0n,(^JXm,, (A5) 

where the notation is: a = {n^;, riy, n^}. For each combination of quantum numbers 
{nx,ny,nz}, the spin part is chosen in such a way that the basis state is an eigenfunc- 
tion of the x-simplex operator Sx = Pe"*'^"'"', where P denotes the parity operator. The 
positive and negative x-simplex operator eigenstates: 

\n,nyn,- i = +) = \nxnyn,) ^ fit) " (-1)"' U)] ' (A6) 
\nxnyn,;t = -) = \nxnyn,) (-l)"^+««+il!l [|t) + (-1)"^ , (A7) 

are related by the time-reversal operator (T = iayKo) 

In^nyTi^] i = -) =f Irixnyn^] i = +) . (A8) 
15 



For the Dirac spinor with positive simplex eigenvalue, the large component corresponds to 
positive, and the small component to negative eigenvalues 



V',(r,+)= \ '''' M . (A9) 




The large and small component are expanded in terms of the basis states Eqs. (IA6I) and 
Mr- +)=Y1 ftM^; +) and g,{r; ") = ^ ^fM^'^ -) ■ (AlO) 

a a 

Positive simplex eigenstates are denoted by and negative simplex eigenstates by \a). 
If the basis states are arranged as: {ai, . . . , ajv/, ai, • • • ajv/}, the x-simplex operator has a 
simple block-diagonal form, whereas the time-reversal operator is skew diagonal: 



and f = . (All) 





If the dimension of each simplex block {i = ±) is M, the dimension of the entire configuration 
space equals 2M. In the present implementation of the model parity is also conserved, and 
this allows a further reduction of the basis to four simplex-parity blocks. For a given maximal 
number of oscillator shells Nmax, the dimension of the HO basis can be determined as follows. 
The states InxnyUz) within a major oscillator shell N are arranged as 





Uy 


nz 








N 





1 


N-1 











1 A^- 1 

1 A^- 1 
A^ 

and the number of 3D HO basis states in the shell A^ then reads 

rist^UN) = {N + 1) + N+{N-1) + --- + 1 = l(iV + l)(A^ + 2) . (A12) 
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Because parity is conserved, the basis can be separated into positive and negative parity 
blocks. The dimension of the each block is determined by summing up the number of states 
in even/odd- TV shells: 



-j^ "'max -j^ 

npo. = -J2i^k + l){2k + 2) = - {ki,, + 1) (/c+,, + 2) (4fc+,, + 3) , (A13) 
^neg. = - J] (2/c + 2) (2 A; + 3) = - (A;_, + l) (/c",, + 2) (4A;-„, + 9) , (A14) 

k=0 

where k^^x — Wmaxf^] and k^^^ = [{N^ax — l)/2], and the square brackets denote integer 
division. 

APPENDIX B: TRANSFORMATION OF THE PRODUCT OF ID HO WAVE 
FUNCTIONS TO THE CENTER-OF-MASS FRAME 

By multiplying the generating function for the Hermite polynomials 

oo ^ 

g{x,p, b) = e^^^/^'-P' = V -p-H^{x/b), (Bl) 

^-^ nl 

n=0 

with the factor obtain the generating function for the HO wave func- 

Vb 

tions: 



_L^-l/.^-xV2,H2x,/t-,^ ^ J2vn{p)Mx,b), (B2) 

where rjnip) denotes 



/ 2" 

^n(p)=p"V^. (B3) 

V nl 

For the product of two generating functions 

g{xi,pi,b)g{x2,P2,b) = l7r"i/2g-5i^(x?+xi)+f 

a new set of coordinates is introduced: 

with an analogous relation for the variables pi and p2- The exponent in Eq. (IB4p can now 
easily be expressed in terms of the new coordinates: 

g{x,,Pub)g{x2,P2.b) = ^^n-'/'e-M^'^^')H{^P+^P)-ip'+P') = g{X , P,b)g{x,p,b). (B6) 
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By using the definition of tfie generating functions 

CXj oo 

g{X,P,b)g{x,p,b) = ^riN{P)(l)N{X ,b)^rinip)(j)n{x,b) 



N=0 



n=0 



the coefficients rjN^P) and rjnip) can be expressed in terms of pi and p2'- 



MP) = 




{P2 + Pi 




M=0 



M 



nl \ nl V nl ^-^ 

m=0 

The product of generating functions Eq. flB7p then reads 

00 I N n I 

n,Af=0 ■ ■ Af=Om=0 \ 



n 
m 



pTp"- 



N 
M 



n 
m 



N-M+m M +n-m 
Pi P2 



With the introduction of the auxihary indices: 

ni = N — M + m and n2 = M + n — m, 
the product Eq. flBlOp can be rewritten in the form: 



g{X,P,b)g{x,p,b)= pTp7 E MXMUi.h)J — 



ni,n2=0 n,N=0 
N n 



N\n\ 



■X 



M=Om=0 \ ^-^7 

One of the Kronecker symbols can be used to eliminate the sum over M: 



g[X,P,b)g[x,p,b)= p^pI- ^ MX MK{xM\ -—5^, 



ni ,n2=0 



m=0 



n,N=0 



I- 

A comparison with the equivalent relation for the product (cf. Eq. ( 1B4[) 





f N \ 










II 






. N — ni + m . 




[ml 



g{xi,pi,b)g{x2,P2,b) = ^ Pi^P2' 



ni,n2=0 



2ni+n2 

ni\n2\ 



n2 

ix2,b) 
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leads to the final expression for the transformation of the product of HO wave functions: 

0„,(xi,6)</.„,(x2,6) =J2M:^n,MX,b)M3^,b) , (B15) 



N,n 

where M^^^ are the coefficients of the transformation 





f N \ 




f „\ 






II 






[ N — rii + m . 




[ml 



^nZ. = ^^^^<.n.,n,Nl^{-ir-\ - II I- (B16) 

For the calculation of matrix elements of the pairing interaction, the center-of-mass and 
relative coordinates are used 

1 1 ~ 

X = -{xi + X2) = —j=X and x = Xi — X2 = \p2x . (B17) 

2 V 2 

The HO wave functions are expressed in terms of X and x: 

6) = (/.^(v^X, 6) = ^^^uEn {flXlh) e-^^"'^'" = ^0iv(X, B) , (B18) 
<^„(x, h) = </.„(x/v^, h) = -^^fnHn (x/V2b) c"^'/^''' = v^^x, 6,) , (B19) 

where we have defined the oscillator lengths B = and br = V2b. Finally, the product 

of two HO wave functions expressed in terms of the center-of-mass and relative coordinates 
reads: 

0ni(Xl,6)0„,(x2,6) = J2K^n,MX,B)M^,br). (B20) 

N,n 
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0.5 60" 0.5 60' 




FIG. 1: (Color online) Self-consistent RHB triaxial quadrupole binding energy maps of the even- 
even isotopes ^^^~^^^Sm in the /3 — 7 plane (0 < 7 < 60^). All energies are normalized with respect 
to the binding energy of the absolute minimum (red dot). The contours join points on the surface 
with the same energy (in MeV). 
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FIG. 3: (Color online) Proton (left column) and neutron (right column) pairing energies of 
I34,i36,i38gj^ in the /? — 7 plane (0 < 7 < 60*^). The contours join points on the surface with 
the same energy (in MeV). 
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FIG. 5: Binding energy per nucleon for the sequence of Sm isotopes, calculated with the 3D RHB 
model and compared to data 



18| . In the inset we display the relative differences (in percent): 
|-^RHBZ _ ^3DRHB^y£;3DRHB^ between the binding energies calculated using the 3D RHB and the 
axial (RHBZ) relativistic Hartree-Bogoliubov models. 
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FIG. 6: (Color online) 3D RHB (filled circles) and axial (RHBZ) empty symbols) results for 
the self-consistent ground-state quadrupole deformations (upper panel), and neutron and proton 



pairing energies (lower panel) of even- A Sm isotopes, 
the separable force (squares) and the Gogny DIS force 



n calculations with axial symmetry both 
J] (diamonds) are used. 
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